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Abstract. The scattering properties of quantum particles on fractal potentials at 
different stages of fractal growth are obtained by means of the transfer matrix method. 
This approach can be easily adopted for project assignments in introductory quantum 
mechanics for undergraduates. The reflection coefficients for both the fractal potential 
and the finite periodic potential are calculated and compared. It is shown that the 
reflection coefficient for the fractal has a self-similar structure associated with the 
fractal distribution of the potential. 



PACS numbers: 03.65.Nk, 05.45.Df 



Submitted to: Eur. J. Phys. 



* To whom correspondence should be addressed (jmonsori@fis.upv.es) 



A transfer matrix method for the analysis of fractal quantum potentials 



2 



1. Introduction 

Both quantum mechanics and elementary solid state physics courses illustrate the energy 
band structure in solids through the one-dimensional Kronig-Penney model that consists 
of a periodic configuration of square- well potentials [H |2]. This problem is usually 
solved by matching the boundary conditions of the wavefunctions at the cell boundaries, 
thus requiring the computation of the determinant of a 4 x 4 matrix [F. Recently, 
some less tedious approaches have been proposed which usually can be readily adapted 
to finite periodic potentials |3]. Among these methods, those based on the transfer 
matrix approach which only uses 2x2 matrix operations in a purely algebraic way 
are the most appropriate ones for beginners Moreover, this method allows to 
introduce a numerical method based on a piecewise constant approximation [H] for a 
general potential, the analysis of defects on slightly aperiodic potentials, and even the 
consideration of more complicated potentials. Among the last, fractal potentials is the 
one which we considered here. 

In recent years the study of fractals has attracted much attention because many 
physical phenomena, natural structures and statistical processes can be analyzed and 
described by using a fractal approach [3 IE] • From a mathematical point of view, fractals 
are self-similar structures obtained by performing a basic operation, called generator, 
on a given geometrical object called initiator, and repeating this process on multiple 
levels; in each one of them, an object composed of sub-units of itself is created that 
resembles the structure of the whole object. Mathematically, this property should hold 
in all scales. However, in the real world, there are lower and upper bounds over which 
such self-similar behavior applies. Fractals are becoming a useful tool to be able to 
model diverse physical systems 01111!, have new technological applications [HI [121 • 

In non-relativistic quantum mechanics, fractals have been used to generate new 
solutions of the Schrodinger equation which are continuous but nowhere differentiable 
wave functions and models for the so-called fractal potentials ^j. Fractal 

potentials allow the analysis of quasi-periodic and nearly stochastic potentials using 
the symmetries induced by the self-similar structure of the potential. Here we consider 
the simplest fractal, the (triadic) Cantor set, as a fractal potential for quantum 
scattering [TH] and tunneUing jTHl [Hj ■ 

In this paper, we present a simple transfer matrix method to obtain the scattering 
properties of Cantor set fractal potentials, which can easily be automated by computers. 
The present method makes easier for the comparison with the finite periodic case and 
shows how the refiection coefficient for the fractal case has a self-similar structure 
associated with the fractal distribution of the potential. Moreoever, it can be easily 
implemented in any computer language, e.g., the Mathematica software package, 
accesible to undergraduate students with only a basic programming experience, so that 
it can be adopted for project assignments in computer physics courses. Starting with the 
implementation of the transfer matrix method for a potential barrier, the extension to 
finite periodic potential is a straightforward one; its extension to Cantor set potentials 
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Figure 1. (a) Local scattering with the i-th interface of the piecewise constant 
potential among the values and Vi. (b) Wave propagation through the constant 
potential Vi as used to calculate the propagation matrix. 



can be based on a recursive implementation, involving the possible improvement of 
student programming skills. Furthermore, these projects can introduce the students to 
the analysis of computational complexity of algorithms, since the Cantor set prefractal 
has a large number of potential barriers and its simulation requires an exponential 
number of matrix products. 

This paper is organized as follows: The next section describes the main facts about 
the transfer matrix method for quantum scattering implemented by piecewise constant 
potentials. In Section 3, the reflection coefficient for the tunnelling on both the finite 
periodic potential as well as the Cantor set pre-fractals are determined and compared. 
Finally, the last section is devoted to the conclusions. 



2. The transfer matrix method in quantum scattering 

Let us consider the one-dimensional, steady-state, linear Schrodinger equation 

fc2 fp I J \ 

where ip{x), m and E are the wavefunction, mass and energy of the particle, respectively. 
The constant h is Planck's constant, and V{x) is the quasiperiodic potential which can be 
represented by a piecewise constant function. Figure Ufa) shows the quantum scattering 
at the i-th interface between two successive constant values of the piecewise potential, 
whose position, without loss of generality, has been taken as a; = 0. In this figure, both 
ipf and are forward and backward plane wavefunctions, respectively, on the region 
where the potential value is V^, and ipi = ipt + V^j"- These wave functions are given by 

V.± = 4±e±i^^^ (2) 



where ki = 2 m (E — Vj) is the local particle momentum, and A- are integration 
constants to be determined by applying the standard boundary conditions at the 
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Figure 2. Piecewise constant potential witli N potential wells Vi with dj as the 
corresponding width. Vq and Vat+i are the surrounding constant potential values 
extended to infinity. 



interface. The continuity of the wavefunctions and the derivatives at the boundary 
are given by 

iJi^iix = 0) = ^i{x = 0), At, + A^, = At + A;, 

^P'i-iix = 0) =tljl{x = 0), ki^i 4ti - Ar_, = ki Af - k, AT, 

where the prime denotes differentiation. Eq. Q is a hnear system of equations which 
can be written in matrix notation as 





and yielding 



(4) 



Here on, the matrix -Djl\ -Dj is referred to as the wave scattering matrix. 

After crossing the i-th interface, the plane wave propagates through the constant 
potential Vi until it finds the next interface at a distance di. Using the notation shown 
in Figure ^b), this wavefunction is given by 

^± = Af e±i e±i ^ = Af e±i ^, (6) 

and a wave propagation matrix Pi can be defined as 

-,1 di 





;i • (7) 



e 

Both the scattering and propagation matrices can be used to solve the general problem 
of the scattering with a piecewise constant potential with potential wells, as shown 
in Figure 121 The successive application of the scattering and propagation matrices yield 
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Figure 3. Finite periodic (a) and Cantor set quasiperiodic (b) potentials where the 
white and black regions denote the potential values and V, respectively. 



and, in the most general form, 



Both the reflection and transmission coefficients of the scattering of a quantum particle, 
incoming from the left, with the A^-well potential is determined by the coefficients of 
the matrix M, 



i I M21 M22 / I 



(10) 



where no backward moving particle can be found on the right side of the potential, so 
Aj^^-^ = 0. The reflection and the transmission coefficients are given by 

R = - — '^-Kt = ^, and T = " ,' — "."t""' = '"^^ o , (Hi 
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3. Presentation of results 

The simplest fractal potential is the Cantor set, shown in Figure El^b), which can be 
obtained by means of an iterative construction. The ffist step {S = 0) is to take a 
segment of unit length. The next one {S = 1) is to divide the segment in three equal 
parts of length 1/3 and remove the central one. In general, at the stage S, there are 2'^ 
segments of length 3""^ with 2*^ — 1 gaps in between. Stage + 1 is obtained by dividing 
each of these segments into three parts of length 3""^"^ and removing the central ones. 
In Figure Efb), only the four ffist stages are shown for clarity. Note that the S'-th stage 
Cantor set pre-fractal can be interpreted as a quasiperiodic distribution of segments 
which can be obtained by removing some segments in a finite periodic distribution as 
shown in Figure a). This distribution at stage pm has (3*-'^ — 1)/2 + 1 potential barriers 
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of length 3~^^, separated by potential wells of the same length, so the "period" of this 
finite structure is A = 2 ■ 3~^^. 

The scattering problem for both the quasiperiodic, Cantor set, pre-fractal potential, 
and the finite periodic potential can be easily solved by means of the matrix transfer 
theory presented in Sec. 2. It is standard to normalize both the energy and the height 
of the potential barrier by the period A, introducing the non-dimensional variables 



(p = A , and cpv = A . 

h n 

Figures m and El show the reflection coefficient, R, for the finite periodic potential 
and Cantor set fractal potential, respectively, around the interval which contains the first 
band gap of the infinite periodic one. Using the standard Kronig- Penney model JHI, this 
band gap can be numerically calculated yielding 3.2519 < (j) < 3.6222 for the potential 
(py = 2. In this energy interval, a Bloch wavefunction does not propagate in a infinite 
periodic potential and, therefore, the transmission coefficient should vanish {R = 1). 
Only evanescent wavefunctions characterized by a complex wavevector, k, are solutions 
of the Schrodinger equation. For this reason, when the number of periods is finite, the 
quantum particle may pass through the potential distribution by the tunnelling effect. 
Figure m shows that the reflection coefficient approaches unity as the number of periods 
in the spatial interval increases, illustrating the process of appearance of the band gap 
of the (full) periodic structure. Although, at the graphical resolution of Figure IH^c) the 
value i? = 1 is apparently reached, the reflection coeficient is always smoUer than unity 
in finity peridodic structures. 

Figure El shows the refiection coefficient for the Cantor set pre-fractal potential for 
S = 2 (top), S = 3 (middle) and 5" = 4 (bottom). It is shown that the refiection at each 
higher stage is a modulated version of that associated with the previous stage. That is, 
the refiection spectrum exhibits a characteristic fractal profile that reproduces the self- 
similarity of the potential distribution. In fact, any wide peak at stage 5* is transformed 
into three narrower and taller peaks at stage S* + 1. Zero refiection from these fractal 
quantum potential occurs at specific discrete energies, while near total refiection is 
possible at other discrete energies. Comparing Figures |31 and an increasing number 
of zeros inside the band gap is observed. These zeros represent resonances due to the 
presence of "defects" in the quasiperiodic potential obtained by removing some segments 
in the finite periodic sequence. 



4. Conclusions 



The transfer matrix method is becoming the standard method for the calculation of 
the tunnelling of quantum particles on constant piecewise potentials because it can 
be used for simple, textbook-like problems and as a numerical method for computer 
simulations. This procedure has been applied to Cantor set fractal potentials, which are 
constant value potentials with support on a Cantor set. For pre-fractals, the S'-th stage 
fractal, the reflection coefficient was numerically calculated and compared with that of 




Figure 4. Scattering reflection coefRcient for the finite periodic potentials of stages 
P2 (a), P3 (b), and pi (c) as a function of the normalized energy for the potential 
(/)v = 2. 




Figure 5. Scattering reflection coefficient for tlie Cantor set pre- fractal potentials of 
stages S — 2 (a), 5* = 3 (b), and = 4 (c) as a function of the normalized energy (p 
for the potential (j>v — 2. 
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a finite periodic potential of tlie same period. Tlie appearance of the first band gap 
of the Kronig-Penney model in the finite periodic potential has been illustrated. The 
reflection coefficient for the Cantor set potential is self-similar. 

The transfer matrix method presented in this paper can be easily adopted in 
computer laboratories for undergraduate quantum mechanics courses, providing a 
powerful method for developing students skill on physics by means of computational 
tools. Furthermore, fractal geometry is a highly motivating topic for the students 
providing a great opportunity to undertake projects closely related to research ones. 
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